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We consider exciton optical absorption in quasiperiodic lattices, focusing our attention on the 
Fibonacci case as typical example. The absorption spectrum is evaluated by solving numerically 
the equation of motion of the Frenkel-exciton problem on the lattice, in which on-site energies 
take on two values according to the Fibonacci sequence. We find that the quasiperiodic order 
causes the occurrence of well-defined characteristic features in the absorption spectra. We also 
develop an analytical method that relates satellite lines with the Fourier pattern of the lattice. Our 
predictions can be used to determine experimentally the long-range quasiperiodic order from optical 
measurements. 

PACS numbers: 71.35+z; 36.20.Kd; 61.44. +p 
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I. INTRODUCTION 

Linear lattices described by means of the binary Fi- 
bonacci sequence have been regarded as one-dimensional 
(ID) quasicrystals and, consequently, they have been the 
subject of intensive theoretical studies. These systems 
lack translational symmetry but, unlike disordered lat- 
tices, display long-range order. The corresponding elec- 
tronic states present rather unusual properties like fractal 
energy spectra and self-similar wave functions. From the 
experimental viewpoint, several electronic properties of 
solids can be inferred from optical characterization tech- 
niques like optical absorption, photoluminescence and 
fluorescence after pulse excitation. Then, a complete un- 
derstanding of the interplay between the electronic prop- 
erties and the underlying quasiperiodic order requires a 
detailed analysis of optical processes. This reason has 
motivated various works dealing with optical properties 
of Fibonacci systems, mainly devoted to Fibonacci semi- 
conductor superlattices 

In previous works we have focused our attention on 
excitons in Fibonacci lattices J^||. In particular, we 
have numerically evaluated optical absorption spectra 
due to Frenkel excitons in self-similar aperiodic lattices 
(Fibonacci and Thue-Morse) pi . We found several char- 
acteristic lines specific of each aperiodic system that are 
not present in periodic or random Frenkel lattices, thus 
being an adequate way for determining the particular 
structural order of the system from experiments. When 
the two kind of optical centers composing the aperiodic 
lattice are very different, the origin of all the absorp- 
tion lines can be successfully assigned by considering the 
aperiodic lattices as composed of two-center blocks, as 
much as in the same spirit of the renormalization group 
concepts developed in the last few years This ap- 

proach neglects the long-range order present in the lattice 
and only takes into account the short-range interaction 
between nearest-neighbor optical centers. These results 
lead, in a natural way, to the question whether long-range 
quasiperiodic order can also be characterized by means 
of optical absorption spectra. It is clear that we must 



search for specific lines due to the long-range order in 
the opposite limit, namely when the two kind of optical 
centers are of similar nature. We will show below that 
this is indeed the case. 

In this work we investigate optical absorption due to 
Frenkel excitons in ID binary systems arranged accord- 
ing to the Fibonacci sequence. We make use of a gen- 
eral treatment which allows us to study the dynamics of 
Frenkel excitons in these lattices, solve the microscopic 
equations of motion and find the optical absorption spec- 
trum. The main conclusions of this paper are twofold. 
First, we show that Frenkel excitons in Fibonacci lat- 
tices with similar optical centers lead to absorption lines 
specific of this kind of ordering and, consequently, not 
present in periodic or disordered lattices. Second, and 
most important, by means of an analytical approach we 
are able to explain the origin of these characteristics lines, 
which are related to the Fourier transform of the lattice. 
Thus we successfully relate an optical property (optical 
spectrum) with an structural one (topology of the lat- 
tice). 



II. PHYSICAL MODEL AND THEORY 

We consider a system of N optically active, two-level 
centers, occupying positions r„ on a ID regular lat- 
tice with spacing unity. Therefore, the effective Frenkel 
Hamiltonian describing this system can be written (we 
use units such that h = 1) 



H = ^ v naia n + ^2 Jnl<4i a l- 

n l^n 



(1) 



Here a^ n and a n creates and annihilates an electronic ex- 
citation of energy V n at site n, respectively. J n \ (n ^ I) is 
the intersite interaction of dipole origin between centers 
n and I, assumed to be of the form J n \ = — J\r n — r/|~ 3 , 
where J is the coupling between nearest-neighbor cen- 
ters. Since J n i is a rapidly decreasing function of the 
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distance between centers, we will omit interactions be- 
yond nearest-neighbors. 

In this work we will be concerned with the Frenkel 
exciton dynamics in systems presenting a quasiperiodic 
distribution of optical centers. The Fibonacci lattice is 
the archetypal example of deterministic and quasiperi- 
odically ordered structure. Any arbitrary Fibonacci sys- 
tem presents two kind of building blocks. In our case, 
we choose those blocks as individual two-level centers 
with on-site energies Va and Vb- The Fibonacci ar- 
rangement can be generated by the substitution rule 
A — ► AB, B — > A. In this way, finite and self-similar 
quasiperiodic lattices are obtained by n successive ap- 
plications of the substitution rule. The nth generation 
lattice will have N — F n elements, where F n denotes a 
Fibonacci number. Such numbers are generated from 
the recurrence law F n — F n —\ + F" n _2 starting with 
Fq = F\ = 1; as n increases the ratio F n -i/F n converges 
toward t= (VE-l)/2 = 0.618..., an irrational number 
which is known as the inverse golden mean. Therefore, 
lattice sites are arranged according to the Fibonacci se- 
quence VaVbVaVaVbVaVbVa ■ ■ ■ where the fraction 
of .B-centers is c ~ 1 — r. 

Having presented our model we now briefly describe 
the method we have used to calculate the absorption 
spectra. The line shape 1(E) of an optical-absorption 
process in which a single exciton is created in a lattice 
with N sites can be obtained as follows ||. Let us con- 
sider the total dipole moment operator T> — X^„(a„+tin), 
where the dipole moment of each center is taken to be 
unity. Here we are restricting to the case of systems 
whose length is much smaller than the optical wave- 
length. Denoting by and Ek the eigenvectors and 
eigenvalues of the Hamiltonian Ji, respectively, the ab- 
sorption line-shape is given by 



7 ^) = ^E l(fc|^|vac)| 2 ^-^), 



(2) 



where |vac) denotes the exciton vacuum. In practice one 
considers a broadened ^-function to take into account 
the instrumental resolution function. Therefore we will 
replace the 5-function by a Lorentzian function of half- 
width a, hereafter denoted by S a . 

A reliable method for determining 1(E) numerically 
involves the consideration of the correlation functions 



Gn(t) = 2J (vac|a n (i)a[|vac), 



(3) 



where a n (t) = exp(iTlt)a n exp(-iTlt) is the annihilation 
operator in the Heisenberg representation. The function 
G n (t) obeys the equation of motion 



dt 



G n (t)=J2 H nlGl(t), 



(4) 



diagonal elements are simply given by — J. The micro- 
scopic equation of motion is a discrete Schrodinger-like 
equation on a lattice and standard numerical techniques 
may be applied to obtain the solution. Once these equa- 
tions of motion are solved, the line shape is found from 
the following expression 



1(E) = - JL jf dt e~ at sin(M)Im G n (t)J 



(5) 



where the factor exp(— at) takes into account the broad- 
ening due to the instrumental resolution function of half- 
width a. 
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with the initial condition G n (0) = 1. The diagonal ele- 
ments of the tridiagonal matrix H n i are V n whereas off- 



FIGURE 1 

FIG. 1. Absorption spectra for a Fibonacci lattice with 
TV = F12 — 233 active centers and on-site energies Va = 4.0 
and Vb = 4.8. Solid line indicates numerical result and 
dashed line the analytical approach. 



III. NUMERICAL RESULTS 

We have solved numerically the equation of motion (^) 
using an implicit (Crank-Nicholson) integration scheme 
]l0| . In the remainder of the paper, energy will be mea- 
sured in units of J whereas time will be expressed in units 
of J -1 . The energy and time scales can deduced from the 
experiment since the exciton bandwidth is 4 J. Fibonacci 
lattices are generated using the inflation rules discussed 
above. In order to minimize end effects, spatial periodic 
boundary conditions are introduced in all cases. We have 
checked that the position and the strength of all lines of 
the spectra are independent of the system size. Hereafter 
we will take N = F\i = 233 as a representative value. We 
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have set Va — 4.0, Vb — 4.8 and J = 1 as typical val- 
ues corresponding to centers with similar characteristics. 
The width of the instrumental resolution was a = 0.15. 
The maximum integration time and the integration time 
step were 50 and 10 -3 , respectively. Smaller time steps 
lead essentially to the same results. 

In the pure A lattice the spectrum is a single 
Lorentzian line centered at E = Va — 2J = 2.0, with 
our choice of parameters. When f?-centers are intro- 
duced according to the Fibonacci rules, a shift of the 
position of the main line towards higher energies is ob- 
served. This shift is also observed in the case of random 
lattices Ql2|. The main line is now located around 
E = 2.30, as shown in Fig. [l]. Besides this main absorp- 
tion line, two satellite lines can be also observed in the 
high-energy part of the spectrum, centered at energies 
= 5.92 and E^ = 4.03. According to the renor- 
malization group techniques developed in Ref. none 
of these satellites can be explained within the so called 
two-center approach. Thus we are led to the conclusion 
that they are due to long-range effects of the quasiperi- 
odic ordering of the lattice. This will be more evident 
after performing the analytical treatment. 



IV. ANALYTICAL RESULTS 

In this section we develop a theoretical approach to 
explain the shift of the main line and the occurrence of 
well-defined satellite lines in the high-energy region when 
.B-centers are introduced quasiperiodically in the lattice. 
To this end we rewrite the system Hamiltonian TL = H.A J r 
Hab, where 

= Va a n a n ~ (alftn+i + a„ +1 a n ), (6a) 

n n 

H AB = ^ (V„ - Va)4«»- (6b) 

n 

Here Ha is the Hamiltonian corresponding to the pure A 
lattice. 

As a first step we consider the pure A lattice whose 
Hamiltonian is TIa- This Hamiltonian with periodic 
boundary conditions can be exactly diagonalized, yield- 
ing the eigenvectors 



W = (jf) E exp ( 2ni ~jPj a ™' vac 



(7a) 



with k — 0, . . . , N — 1. The corresponding eigenvalues 
read 



E k = V A - 2 J cos 2tt 



N 



(7b) 



Therefore, the dipole moment matrix elements within 
this approximation satisfy the relation (fc|Z>|vac) = Ska- 
Consequently, one obtains that 1(E) = (l/N)8 a {E — E ). 



This result shows that the pure A lattice presents a single 
Lorentzian absorption line centered at E = Va — 2 J, as 
we have mentioned above. 

To develop an approach able to explain the results in 
Fibonacci systems one must explicitly consider the ef- 
fects of the topology through the term Hab- To carry 
out such an approximation we have used the nondegen- 
erate perturbation theory. Let \k^ AB ^) the perturbed 
eigenvector describing the exciton state when the influ- 
ence of the term Hab is taken into account. To evaluate 
the perturbed dipole moment matrix element, we expand 
\kS AB ^) in the basis of the unperturbed eigenvectors \k) 
given in ( |7a| ) 



i^ B >> = i*>+E 



(i\n AB \k) 

Eh — Ei 



\l). 



(8a) 



whereas the perturbed eigenvalues are given by 



E[ AB) =E k + (k\H AB \k) = E {AB) + 4 J sin 2 ( n- j , 



where E. 



(AB) 




(8b) 

cV B + (1 - c)V A - 2 J. Notice that Hab 
shifts all levels by the same amount. This shift is given 
by (k\T-L A B\k) = c(V B - Va) for all k. With the chosen 
parameters it amounts 0.30, in perfect agreement with 
our numerical results shown in Fig. [l]. 
Inserting (||) in (Q) one finally obtains 

1 1 N ~ 1 

1(E) = -8 a (E - 4 AB) ) + \F(k)\ 2 S a (E - E {AB) ), 



fc=i 



where for brevity we have defined 



F(k) 



(o\n AB \k) 

4Jsin 2 (7rfc/iV)' 



(9) 



(10) 



Therefore, the spectrum consists of a main Lorentzian 

line centered at E [AB) = cV B + (1 - c)V A - 2 J = 2.30, 
in agreement with the virtual crystal approximation, and 
several satellite lines at E^ AB ^ . The most remarkable fact 
of expression (|^), however, is that it relates the finer de- 
tails of the absorption spectrum, an experimentally mea- 
surable magnitude, with the Fourier transform of the lat- 
tice describing the topological long-range ordering of op- 
tical centers. To demonstrate this point, we recall that 
unperturbed eigenvectors are orthogonal. Thus 

(0|Wab|*) = (0| ( F « " VA)4«n|fc) 
n 

n ^ ' 

which is nothing but the aforementioned Fourier trans- 
form. A comparison of the perturbative prediction (||) 
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with the numerical results is given in Fig. [j], where an 
excellent agreement is achieved in the weak-perturbation 
limit we have considered. In particular, our analytical 
approach predicts the existence of two satellite lines at 
EM = 5.78 and = 4.13 

In order to get a deeper insight into the problem, 
we relate our results with those reported in the vari- 
ous works dealing with the Fourier transform of the Fi- 
bonacci sequence, some of them in connection with X- 
ray diffraction in aperiodic semiconductor superlattices 
p3| , p^| . The Fibonacci lattice is known to exhibit a 
Fourier spectrum displaying well-defined Bragg peaks. 
The Fourier intensity |(0|7Yab|^)| 2 consists of a series 
of peaks located at values of the momentum 2ir(k/N) of 
the form r p , where p is an arbitrary integer. Each peak 
of the Fourier intensity leads to a large contribution to 
the sum (^|) and, consequently, to a well-defined satellite 
line in the absorption spectrum. 
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FIGURE 2 

FIG. 2. Plot of |F(fc)| 2 as a function of the ratio k/N for Fi- 
bonacci lattices, showing the ocurrence of well-defined Bragg 
peaks, labelled according to successive powers of r. 

Figure | shows \F{k)\ 2 as a function of the ratio k/N 
evaluated from expression (|l(]), the results being inde- 
pendent of the system size. A sequence of Bragg peaks, 
arranged according to successive powers of r, are clearly 
observed in this plot. The stronger one, which is re- 
sponsible of the highest-energy satellite line centered at 



£?(!) = 5.78, is located at k p /N 



0.38 



(P = 2). 



Therefore the mode is located close to the top of the 
excitonic band (k/N = 0.5). On decreasing the momen- 
tum, the next Bragg peak is located at k p /N ~ 0.23 ~ r 3 
(p = 3), leading to an energy value 4.13, in perfect agree- 
ment with the energy of the second satellite peak in the 
high-energy region of the spectrum observed in Fig. [l]. 
The analytical treatment predicts the existence of other 



satellite lines with smaller energy, but they are hidden 
by the main absorption line. 

Our results are suitable for a direct physical inter- 
pretation. Each satellite line of high energy is caused 
by the coupling of two modes, namely the lowest-lying 
and k p , through the topology of the quasiperiodic lattice. 
Notice that different arrangement of centers would lead 
to a completely different Fourier intensity |(0|7iAs|fc)| 2 - 
Thus, the exciton acts as a probe of the long-range or- 
der of the quasiperiodic lattice. This is one of the main 
results of the present work since it provides us with a 
useful method to be applied in experimental situations. 



V. CONCLUSIONS 

In summary, we have studied numerically the absorp- 
tion spectra corresponding to the Frenkel-exciton Hamil- 
tonian on self-similar quasiperiodic Fibonacci lattices. 
Besides the main line, which is shifted respect to that 
of the pure lattice, we found several satellite lines in the 
high-energy region of the absorption spectra. We have 
also developed an analytical method which successfully 
explains not only the shift of the main line and the posi- 
tion of the satellite lines but also their shape. Our anal- 
ysis clearly indicates that each satellite line is directly 
related to a particular Bragg peak of the Fourier trans- 
form of the underlying lattice. We have realized that each 
satellite line comes from the coupling between the lowest- 
lying mode and a particular mode lying close to the top 
of the excitonic band. This relationship surely should fa- 
cilitate future experimental work on optical properties of 
molecular systems with long-range order. 
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